physics/0302070 
SINP/TNP/03-02 

Reflection and refraction of electromagnetic waves at the plane boundary between two 

chiral media 

Jose F. Nieves 

Laboratory of Theoretical Physics, Department of Physics, P.O. Box 23343 
University of Puerto Rico, Rio Piedras, Puerto Rico 00931-3343 

Palash B. Pal 
Saha Institute of Nuclear Physics, 1/AF Bidhan-Nagar 
Calcutta 700064, India 
(Dated: February 2003) 

This work is concerned with the propagation of electromagnetic waves in isotropic chiral media 
and with the effects produced by a plane boundary between two such media. In analogy with the 
phenomena of reflection and refraction of plane electromagnetic waves in ordinary dielectrics, the 
kinematical and dynamical aspects of these phenomena are studied, such as the intensity of the 
various wave components and the change in the polarization of the wave as it crosses the boundary. 
As a prerequisite of this, we show that the plane wave solution must be written as a suitable 
superposition of the circularly amplitudes on both sides of the interface, we elucidate which is the 
appropriate set of conditions that the solution must satisfy at the boundary, and we set down the 
minimal, and complete, set of equations that must be solved for the coefficient amplitudes in order 
to satisfy the boundary conditions. The equations are solved explicitly for some particular cases 
and configurations (e.g., normal incidence), the salient features of those solutions are analyzed in 
some detail, and the general solution to the equations is given as well. 

I. INTRODUCTION 

The problems related to the propagation of electromagnetic waves in a homogeneous medium, including those 
of reflection and refraction at the interface between two media, form a traditional set in classical electrodynamics, 
and are discussed in any textbook on the subject [|. A less familiar one is the problem of the propagation of an 
electromagnetic wave in a chiral medium. In such a medium, the dielectric constant has a different value depending 
on the polarization of the wave, which implies that the two circularly polarized waves propagate with different speeds. 
A chiral medium, also called "optically active" , can exist for various reasons. For example, a material medium can 
become optically active under the influence of an external field, such as a magnetic field (the well-known Faraday 
effect discussed in many plasma physics textbooks 0]) or a (control) plane electromagnetic wave Other kind of 
media exhibit the phenomenon naturally, such as a gas of chiral molecules Q , or a neutrino gas 0, |g| that contains an 
unequal number of neutrinos and antineutrinos. Recently, man-made composite materials have been produced with 
similar properties 0. 

Various approaches have been proposed in the literature to study the electromagnetic properties of such optically 
active media0, 0, 0, IToL ITU IT^L [la . IT3. ITU [T6l | . As we showed in Ref . ^3 > a particularly simple and general method 
consists in parametrizing the induced charge and current in the most general way consistent with general principles, 
such as the conservation of the electromagnetic current and, if applicable, the isotropy and homogeneity of the medium. 

However, the formalism developed in Ref. | 17| concerns a homogeneous medium only. The problems associated 
with electromagnetic wave propagation in inhomogeneous chiral media were not considered in our earlier work. An 
important example of an inhomogeneous medium is provided by two semi-infinite homogeneous media, separated by 
a plane interface. As is well-known, if an electromagnetic wave is incident on this interface, part of it is reflected back 
into the original medium, and a part is refracted into the other medium. The main object of this paper is to extend 
our minimal and general parametrization to treat this kind of problem. In particular, here we study the problem of 
reflection and refraction at the interface in detail, for the case in which one or both of the semi-infinite media on 
either side of the boundary exhibits the chirality property. 

The paper is organized as follows. In Section [HJ we summarize the parametrization of the induced charge and 
current for a homogeneous chiral media that was introduced in Ref. [lTj . We then extended it to treat inhomogeneous 
media, and we consider in particular the system that consists of two semi-infinite homogeneous media separated 
by a plane boundary, which is what concerns us in the present work. To set up the stage for the study of wave 
propagation across the boundary between two chiral media, in Section UTTl we review briefly the salient features of the 
propagation of a plane wave in a homogeneous chiral medium. Then in Section Hvl we set up the problem of reflection 
and refraction at the plane interface between two chiral media, in terms of a set of equations that must be solved 
for the various coefficient amplitudes that appear in the suitably decomposed plane wave solution of the Maxwell 
equations. In Section El we solve those equations to obtain the amplitudes of the reflected as well as refracted waves. 
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The corresponding solutions for various particular cases (including simpler configurations such as the case of normal 
incidence), and some of their features, are analyzed in some detail before the solution for the general case is given. 
Section IvTI contains some concluding remarks and outlook. 



II. PARAMETRIZATION OF THE INDUCED SOURCES 



Here, and in the rest of the paper we deal exclusively with materials that are isotropic, and we also assume that the 
field strengths are such that the material responds linearly. We use the Heaviside-Lorentz system of units throughout. 



A. Non-dispersive medium 

For the sake of orientation, let us consider first a non-dispersive medium. Our argument, advocated in Ref. |l7j . 
was to use the Maxwell equations themselves as a guide and write the induced charge and current densities in the 
usual form 

And = -V-P, 
dP 

Jind = -jfi + V X M, (2.1) 

but the relationship between the electric and magnetic polarization P and M being more general than the usual one. 
For the example we are considering, the relation is 

P = (e- 1)1? 

M = (1 -fi-^B-CE, (2.2) 

which is consistent with all the general principles stated above. In particular, the isotropy of the medium implies that 
the quantities e, /i, C are scalars rather than tensors, and the non-dispersive nature of the medium implies those same 
quantities are constants independent of the coordinates. As a consequence of the homogeneous pair of the Maxwell 
equations, 

V-B = 0, (2.3) 
3B 

Vx£ = --, (2.4) 
an equivalent parametrization of the induced charge is furnished by writing 

P = (e-l)E + CB, 

M = (1-^)5, (2.5) 

instead of Eq. (|2.2() . In any case, 

Pi„d = -(e-l)V-f, 

jiud = (e - 1)-^- + I 1 - -J V X B - CV X E . (2.6) 

A question that arises naturally is whether the expressions for P and M can both contain two terms, one proportional 
to E and one to B, namely, 

P = (e-l)£ + Ci£, 

M = (1 - ^)B - ( 2 E . (2.7) 

Indeed, this has been proposed in the literature 20| . However, the discussion above makes it clear that such a 
parametrization is redundant, since the expression for the induced sources would again be given by Eq. (53^. with 

C = Ci + C 2 , (2.8) 
and therefore £ is the only physically relevant quantity and not (12 separately. 
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B. Dispersive but homogeneous medium 



While these arguments refer to a non-dispersive medium, similar considerations, and the same conclusions, hold for 
a dispersive, but homogeneous medium as well. In a dispersive medium, the relation between P, M and the fields is 
non-local. Therefore, expressions like eE should interpreted in the form 



e£-> eE 



d 4 x'e(x,x')E(x') 



(2.9) 



and similarly with the other terms involving /i and C (in this section, we will use the four-dimensional notation x = (t, x) 
in the arguments of the functions for compactness.) However, the assumption that the medium is homogeneous implies 
that the functions e, /i, £ depend only on the relative distance x — x 1 and, furthermore, the causality principle requires 
the time dependence to be on the variable t — t' and not on t and t' separately. In this case the relations such as those 
in Eq. (|2.9(l become simple algebraic relations in Fourier space and, as we showed in detail in Ref. [12]], the argument 
and the parametrization we have used above continue to hold, if every symbol is taken to be the Fourier transform of 
the corresponding coordinate space variable, and if we make the identification 



V 

d_ 

at 



ik . 



Thus in Fourier space, the induced sources are parametrized as 

And(^,fc) = -i{e- l)k- E, 

ji n d(w,fc) = — i(e — 1)luE + i(l — /i _1 )fc x B — iQk x E , 
and the corresponding expressions in coordinate space are 



And = -V-[(e-l)E], 
Jind = ^[(e-l)£]+Vx 



1 

1 



D 



- V x ((E) 



It should be noted that the last term in Eq. (|2.13[1 can be written in terms of B using the relation 

d((B) 



V x ((E) 



dt 



(2.10) 

(2.11) 

(2.12) 
(2.13) 

(2.14) 



which as a consequence of Faraday's law, as can be seen very simply in Fourier space. In coordinate space, it is proven 
by noticing that, for a homogeneous medium as we are considering, we can use manipulations such as 



^- I d 4 x((x-x')E(x') 

UXi 



E(x') 



(2.15) 



where we have integrated by parts. 

Thus, in absence of external sources, the Maxwell equations in a homogeneous chiral medium consist of the pair 
given in Eqs. I|2.3|) and l|2.4|) . augmented by 



V x 



V ■ (eE) = , 

1 - » -A d , -, 
B + (E)- —(eE) = 0. 



s /j j dt 

If the external sources are non-zero, they would appear on the right hand sides of these two equations. 



(2.16) 
(2.17) 



C. An inhomogeneous medium 



In an inhomogeneous medium, the function e(x, x') introduced in the previous section, and similarly fi and depend 
on the coordinates x and x' separately, and not just on the relative coordinate x — x'. In this case, the language of 
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the Fourier transforms is not useful. Moreover, whichever way we would like to look at it, the manipulations that led 
us to Eq. 1)2. 14[) no longer hold, and therefore in this case 



V x ((E) ± 



dm 

at 



(2.18) 



The important implication of this for us is that the question that we posed ourselves in Eq. (|2.7[) comes to haunt us 
again. In the present context, it amounts to asking whether instead of Eqs. H2.12fl and (|2.13l) we can write 



Pind = -V- [{e-l)E}- V- (fjB) 



Jind 



[(e-l)£]+V 



1- - ) B 



-Vx (CE) + -(r 1 B). 



(2.19) 
(2.20) 



The issue now is not whether the two parameters £ and rj are redundant and whether one can be absorbed in the 
other. They are not redundant; the last two terms in Eq. (|2.2(J|) are not proportional to each other in general, and 
the second term in Eq. (|2.19() is not zero. Thus, £ and rj have a different and separate physical meaning, and they 
describe different physical effects. 

The question of which physical systems are described by this parametrization is one that we cannot consider here. 
What we can state is that there are no general physical principles that exclude that possibility. Our parametrization 
is both general and minimal, subject to the restriction of linearity and isotropy that we have assumed. However, as 
we have argued, such systems must necessarily be inhomogeneous. As interesting as the exploration of these systems 
would be, further considerations along these lines in a general way is outside our scope. 



D. Two homogeneous chiral media with an interface 

We finally get to consider the kind of system that concern us in this work, namely, two semi-infinite homogeneous 
media, separated by a plane interface. As we argue below, the proper way to parametrize this kind of system is by 
putting rj = in Eqs. I|2.19|) and (|2.2U|) . That is, the Maxwell equations are the same as stated in Eqs. (|2.16|) and 
(|2.17|) . The argument is actually subtle, but simple. 

Let us take the plane boundary to be the plane defined by z = 0. In the idealized limit of taking the plane boundary 
to be infinitesimally thin, then we can represent f] in the form 

fj = r)6(-z)+ri'0{z). (2.21) 

However, on either side of the boundary (that is, for z ^ 0), the function is zero since we are taking both semi-infinite 
media to be homogeneous. Therefore, r) = everywhere. 

The exception to this argument arises if the situation is such that the plane boundary cannot be taken to be 
infinitesimally thin, but instead we must take into account the fact that it has some non-zero thickness 5. Certainly, 
in the transition from one material to the other, the medium is inhomogeneous and in that region fj can be non-zero. 
If we look at the physical effects that occur in the bulk of either media at distances, say \z\ ~ L with L S, then 
the effects due to the non-zero value of 77 in the transition region are of order S/L. In other words, they are surface 
effects. At points that are close to the plane boundary (\z\ ~ 5) such effects could be observable. But deep in the 
bulk of either media the effects become negligible as S/L — > 0, as it should be since in that case the transition region 
can be idealized to be infinitesimally thin and then the representation given in Eq. I|2.21|l is again a valid one. 



III. WAVE PROPAGATION IN AN OPTICALLY ACTIVE HOMOGENEOUS MEDIUM 



A. Dispersion relation 



For illustrative purposes, let us consider a steady-state (monochromatic) electromagnetic wave propagating along 
a given direction, which we take to be the z axis. The vector potential is of the form 

A(z,t) = e~ iujt a(z) (3.1) 

where a(z) must be chosen such that the vector potential A(z, t) satisfies the Maxwell equations with the appropriate 
boundary conditions. Putting this into Eqs. <|2.16|) and 12.17fl . it follows that the wave equation is satisfied only for 



5 



the circularly polarized waves whose polarization vectors are defined by 

e± = -y=(ux ± iuy) , (3.2) 

where we denote by u x , u y and ii z the unit vectors along the x, y and the z directions, respectively. Further, these 
polarization states satisfy the dispersion relations 

=n CT (w,A CT ), (3.3) 

where the refractive index functions n± that appear in Eq. (|3.3|1 are given in terms of the dielectric function (e), 
magnetic permeability (/z) and activity (£) constants of the medium bv[l7l| 

2 - 1 1 i ^ fa A\ 

v ± = — = — ± —F7 ■ ( 3 - 4 
n± t[i eK 

While the index of refraction can in general be complex, even without the C term, in this work we consider only 
systems for which it is real; i.e., the medium is non-absorbing. This implies, in particular, that C is purely imaginary. 
This, in turn, has implications on the microscopic properties of the system with respect to the discrete space-time 
symmetries, which were considered at length in Ref. 



>rop 

m 



B. Linearly polarized wave 

For a plane wave moving in the z-direction that is linearly polarized (i.e., it contains equal admixture of the two 
circular polarizations), the vector potential is of the form shown in Eq. I|3.1|) . where 

a{z) = A(e tzK +e + + e lzK -e_) . (3.5) 

In writing Eq. Q3.5JI we have implicitly chosen the origin and orientation of the coordinate system in such a way that, 
at z — 0, the linear polarization vector of the wave points along u x , or the x-direction. By simple algebra, and using 
Eq. 13.2J1 . A(z, t) can be expressed as 

A(z, t) = Ae- luJt e lA{z) (u x cos 9(z) + u y sin 6{z)) , (3.6) 
v2 

where 

6(z) ee l(K_-K + )z, 

A(z) ee ±(K + + K_)z. (3.7) 

Thus, at a given distance z — d, the polarization vector of the wave points at an angle given by 9(d) relative to the 
x-axis, the phenomenon known as optical rotation. 

However, it should be noted that for a fixed point z, the direction of the polarization vector is fixed and does not 
change with time. This contrasts with what happens for a circularly polarized wave, e.g., 

A(z, t) = A+e- luJt e lzK +e+ . (3.8) 

In this case the direction of the electric field is given by the vector 

e(z, t) = u x cos(ujt — zK+) + u y sin(wf — zK+) , (3-9) 

where, for simplicity, we have taken A + to be real, but similar considerations apply otherwise. For a fixed point z, 
the polarization vector rotates in a circle, in the plane perpendicular to the direction of propagation. 
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C. Elliptically polarized wave 

The obvious generalization of the above special forms is the superposition 

A{z,t) =e- <wt [A + e lzK +e++ A_e tzK -eJ\ , (3.10) 

which contains an unequal mixture of both circular polarizations. A simple geometric representation of this solution 
is given by the following construction. At the point z = 0, the direction of the electric held is given by the vector 

e(z = 0,t) = (A + + A^)u x cosujt + (A + - AJ)u v sinut . (3-11) 

That is, the polarization vector rotates in an ellipse, with the axis of the ellipse lying along the x and y directions. 
For an arbitrary point z ^ 0, we write Eq. i|3.10[l in the equivalent form 

A(z, t) = e -^V A(z) \A + e- ie(z) e + + A_e ie(z) e 1 , (3.12) 

where 6{z) and A(z) are defined in Eq. I|3.7|l . The direction of the electric field is easily obtained from this to be 
given by the vector 

e(z, t) = (A+ + A-)cos(wt - A(z))u x i + (A+ - A-)sin(ut - A(z))u v > , (3.13) 
where we have defined 

u x : = cos 9(z)u x + sm9(z)uy 

u y ' = — sin 8(z)ii x + cos 8(z)u y . (3-14) 

Therefore, at a point z = L the polarization is also elliptical, but in contrast to the standard elliptical polarization 
phenomena in ordinary media, the ellipse itself has rotated by the angle 0(L). In addition, the polarization angle 
acquires a negative phase A(L) relative to the x' axis. 

In the next sections, we extend these examples to the case in which the wave crosses the plane interface between 
two media, each of which is individually described as we have explained above. The treatment of that problem must 
necessarily take into account the kinematical and dynamical aspects of the reflection and refraction of the wave at 
the plane boundary, among other issues. 

IV. REFLECTION AND REFRACTION OF A CIRCULARLY POLARIZED WAVE 

A. Kinematics 

We consider an electromagnetic wave incident on a plane boundary as shown schematically in Fig. 2] The homo- 
geneity of the space in the xy directions allows us to seek the (plane wave) solutions in the form 

e-*"V Sj - ,s \£(jb). (4.1) 
For ordinary (non-chiral) media, a suitable ansatz is 

A(z) = \de lk ^ + be~ lk « z ] e(-z) + \ce zk '« z ] Q{z) . (4.2) 



The form chosen in Eq. H4.2J) corresponds to the physical situation in which the wave is incident from the left (with 
no reflection from the far right-hand side, as it should be in a semi-infinite medium), as depicted in Fig. ^ 

But here we encounter the first departure from the standard treatment when we consider chiral media. Namely, 
while we can set the incident component to consist of only one of the two polarizations, the reflected and the refracted 
ones will in general consist of a superposition of the two propagating modes, which in the present case have different 
wavelengths. For example, if we decompose the refracted wave vector into its perpendicular and parallel components 
(relative to the z = plane) k' — (k±,k'u), then for fixed values of u> and k± there is actually a different value of 

the parallel component for each polarization, given by a/ K' t 2 — k\ , so that the two polarization components travel in 
different directions. 



(I) 



2 = 



(II) 




FIG. 1: Schematic diagram showing a plane electromagnetic wave with wave vector k incident on a plane boundary at z — 0, 
giving rise to a refracted and a reflected wave with wave vectors k' and k", respectively. 



With this in mind, we therefore write our proposed solution corresponding to a wave with definite polarization 
er = ± incident from the left, in the form 



A(x,t) = e~ iuJt AiQ(-z) + A u Q(z) 
where Aj, which contains the incident and the reflected components, is taken to be 



r=± 



while the refracted component An is 



.4 



ii 



(4.3) 



(4.4) 



(4.5) 



T = ± 



The wave vectors in the expressions for these waves can be written as 

% = (kx,K\\), 
k T = (fcj_ , k T || ) 



(4.6) 



sharing a common value of the perpendicular component. For the same fixed value of u>, we denote by K a {uj) and 
K' a (oj) the solutions of the dispersion relations in the two regions characterized by refractive index functions n± and 
n' ± . Making reference to the angles shown in Fig. ^ the parallel components for the different wave vectors can be 
written as 



K a cos a 
K' T cos a' T , 
K T cos a" , 



(4.7) 



where 



kL 



K, 



l 2 



k 2 



K 2 T -k\ 



(4.8) 
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In order to specify the polarization vectors uniquely, we define e-y to be in the direction of k± x u z . Making reference 
to Fig. ^ it points perpendicular to the page, away from the viewer. We take that to be also the direction of the x 
axis, and define the y axis by 

u y = u z x u x , (4.9) 
which is the vertical direction in the representation of Fig. ^ Then, for the polarization vectors, we take 

e[ = e'l = ei = u x , (4-10) 

while 

e 2 = k X ei , 
e 2r = k T x e\ , 

e 2 ' T = K x ei , (4.11) 

where fc, k T and fc" are the unit vectors along k a , k' T and k", respectively. Notice that, while the direction of the 
incident wave is fixed (by us), the reflected and refracted components consist each of two polarized waves that travel 
in different directions, as we have already remarked. With the choice of axis specified above, 

k = cos au z + sin a u y , 



k T = cos a T u z + sin a' T u y , 

k" = — cos a" u z + sin a" u y , (4-12) 



and whence, using Eq. I|4.11|l . 



e 2 = — sin au z + cos a u y , 
e' 2r — — sin a' T u z + cos a' T u y . 

e 2r = — sin a T u z — cos a" u y . (4-13) 

The circular polarization vectors that enter in Eq. (|4.3(l are defined by 

1 . . x 

e CT = y|( e i + i0_ e 2 ) , 

e' T = -^=(ei +ire 2T ), 

= ^(e 1+ *re 2 ' T ), (4.14) 

and they satisfy 

^•fe^ = 0, 
e' •£ = 0. 



Together with Eq. I|3.3fl . in this way it is ensured that the functions 



e 



-iwt 



A x (X = 1,11), (4.15) 



satisfy the Maxwell equations in each region. 

The fact that the wave vectors have the same transverse component, as indicated in Eq. I|4.6[l . implies the familiar 
relationship (SneH's law) between the angles of incidence and refraction, 



sin a' K a n a 



sin a K' T n' T 
and the analogous relation for the reflected wave 

a T = a (t = a 
sin a" K a n a 



(4.16) 



sin a K T n T 



(r^a). (4.17) 
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Thus, for a given value of u> and a, so that the wave vector fe CT of the incident component is completely specified, 
the wave vectors fc" and k' T of the reflected and refracted components are also completely determined, as well as the 
polarization vectors. The only quantities yet to be determined are the amplitudes a' T , a" in Eq. Ij4.3|) . 

From Eqs. (|4.16|l and (|4.17(l the familiar effects such as total internal reflection can be deduced, in the form that 
they apply to the present situation. Thus, for example, if the dielectric constants are such that K a > K' T , there is a 
maximum incident angle a Tma x, determined by 




(4.18) 



above which the refracted wave with polarization r propagates parallel to the surface (k' T u = 0) but does not penetrate 
the region II. 

In general, the picture that emerges is the following. The incident wave, which has a definite circular polarization, 
is split at the boundary into a reflected and refracted wave, each of which is a superposition of the two circular 
polarizations. However, the two circularly polarized modes do not travel in the same direction as a consequence of the 
fact that their wave vectors have the same transverse component but a different longitudinal one. The determination 
of the relative admixture of the two polarization modes in the reflected and the refracted wave is precisely one of the 
dynamical questions that we will address below. 



B. Dynamics 

Some of the dynamical issues that we wish to consider involve (i) finding the amplitudes of the reflected and 
refracted components of the wave for a given incident amplitude; (ii) determining how the polarization of the wave 
changes as it crosses the boundary. In contrast with the kinematical properties that we considered in the previous 
section, these dynamical aspects depend on the specific nature of the boundary conditions at the interface. The latter 
are determined by the requirement that the solution given in Eq. (|4.3I) satisfies the Maxwell equations at the interface. 



1. Boundary Conditions 

As we have already mentioned, the dispersion relations given in Eq. I|3.3[) together with the conditions given in Eqs. 
I|4.15(l imply that the function A(x, t) defined in Eq. (|4.3(1 satisfies the Maxwell equations in each region separately. 
However, we have yet to ensure that the equations are satisfied at the boundary between the two regions itself. This 
requirement yields further conditions. Are they the same as for the ordinary equations? Or does the £ activity 
constant term modify them in any way? Answering this question is our next hurdle. 

To begin, we calculate the electric and magnetic field associated with the vector potential given in Eq. Ij4.3|l . The 
standard formulas 



E 



OA 

B = V x A. 



yield 



E = e 



B 



E I @(-z) + E II e(z] 



BM-z) + B u Q(z) + 8(z)u z x {A n - Aj) 



where 



Ex 
B x 



iujA x , 
V x A x 



(X = I,II) 



(4.19) 

(4.20) 
(4.21) 

(4.22) 



Since the fields must be finite everywhere (otherwise the Maxwell equations would be not be satisfied) , the last term 
in Eq. (|4.21() must be zero, which yields the condition 



(«re T 







(4.23) 
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on the polarization vectors. 

The remaining relations follow from the Maxwell equations themselves which were given in Eqs. (I2.3|) . (|2.4I) . (|2.1()|) 
and H2.17|l in the Introduction. The argument to obtain the implied conditions is similar to the one used to arrive at 
Eq. (I4.23fl . For example, consider Eq. I|2.16|l . Writing 



IE = e 



we then calculate 



eEiQ(-z) + e'E u Q{z) 



9(-z)V ■ {eEi) + 0(z)V • {e'En) 
+ 5{z)u z ■ (e'En - eEj 



(4.24) 



(4.25) 



The first two terms on the right-hand side are identically zero since there is no free charge in either region. Thus, the 
condition that the divergence of eE is zero everywhere requires that the coefficient of S(z) in Eq. (|4.25|l be zero. This 
yields the relation 



(ea"e" 



= 0, 



(4.26) 



where we have used the expressions for Ei, Ejj that follow from Eqs. I|4.3() . 14.4fl . (|4.5I) and (|4.22() . 

The conditions that follow from the remaining Maxwell equations are obtained similarly. In general, by means of Eqs. 
(|3.3|) and the definition of the corresponding polarization vectors given above, it is ensured that the proposed solution 
satisfies the equations in each region separately. Thus, when the solution is substituted in the Maxwell equations, only 
the terms that arise from the derivatives acting on the Theta function, which are proportional to the delta function, 
are not automatically zero. Demanding that they be zero as well yields the required relations. Proceeding in this way 
we find that Eq. I|2.4|) yields Eq. (|4.23|) again, and the remaining two conditions that complement Eqs. (|4.23|) and 
are 



u z x 



aa a e a K a + ^ r (K r a'X - K' T a' T e' T ) 

T 

aa^evaea + r (a"ew T e" — a' T e'v' T e' T ) 



= 0. 



= 0. 



In writing these conditions, we have used the relation 

%k(j X 0(j — — (T&o 



(4.27) 
(4.28) 

(4.29) 



and the definition of v a from Eq. (|3.4|l , as well as the analogous ones for the corresponding primed and doubly-primed 
quantities. In addition we have used the fact that, for any wave (incident, refracted and reflected), the product of its 
wave number times its velocity velocity equals u), which is the same for all of them. 



V. SOLUTION TO THE BOUNDARY CONDITIONS 



A. The independent conditions 

Our task is to use the conditions given in Eqs. I|4.23[) . I|4.26|l . I|4.27[) and (|4.28() to obtain, for a given incident 
amplitude a CT , the reflected and transmitted amplitudes a' T ,a". There are thus a total of four unknown variables to 
be solved for. On the other hand, at first sight, the number of conditions seems to be six. While Eqs. I|4.26|l and 
(|4.27|l are scalar equations, the relations in Eqs. 14.23fl and (|4.28|l are vectorial. Being a two-dimensional problem, as 
we have seen, each of the vector relations yields two conditions, for a total of six. Therefore, the solution for a' T , a" is 
algebraically overdetermined, and a non-trivial solution exists only if additional auxiliary conditions are satisfied. As 
we will see, these are just the Snell-like relations given in Eqs. I|4.16|l and (|4.17|) . and the dispersion relations of Eq. 
(|3.3|) . Put in another way, when Eq. ll-'S./ifl as well as Eqs. (|4.16|) and (|4.17|) are satisfied, two of the six relations implied 
by the set of boundary conditions are redundant, leaving just four independent equations for the four unknowns a' T , a". 
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Before going to the details of the problem, it is useful to recollect the corresponding situation in the case of ordinary, 
non-chiral media. In that case, for a given value of the transverse component of the wave vector and a given frequency, 
the longitudinal component of the wave vector in each of the reflected and refracted waves has a unique value. Thus, 
the solution written down in Eq. (|4.3|l indeed collapses to the form given in Eq. 1)4. 2|) . Nevertheless, the problem 
is similar to the one stated above. Namely, for a given incident amplitude a, we wish to know what are the the 
coefficients b and c. Being a two-dimensional problem as already argued, then for an arbitrary (but definite) a in the 
ei.2 plane, there are four unknowns, represented by the two components of b in the e" 2 plane and the two components 
of c in the e'-y 2 plane. Further, the problem being a linear one (the boundary conditions are linear equations in the 
amplitudes), we can proceed to find the particular solutions corresponding to the case a — e.\ and, separately, a — e 2 , 
and the general solution for an arbitrary choice of a is obtained by linear superposition of those particular ones. This 
is actually the procedure followed in Jackson's book [l|. The fundamental reason why a solution can be obtained at 
all, is the fact that there are four independent algebraic relations (that follow from the boundary conditions) among 
the four unknowns the we have identified. 

Returning to our problem, these same general principles hold as well. Namely, while some features and details of the 
solutions are not applicable or relevant, the linearity property, superposition principle and the fact that the problem 
is algebraically well defined, also apply. The fundamental difference, as far as the algebraic details and manipulations 
is concerned, is that in the present case it is not useful to build up the general solution by superposing the particular 
solutions obtained for the cases a — and decomposing the reflected and refracted waves in terms of the linear 
polarization components. In our case, the circular polarization basis is, as we have shown, the appropriate one to use. 

For the remaining vector algebra we use the relations 

» » 1 

u z ■ e a = -r=icrsma, 

. 1 . . , 
u z ■ e T = -=tTsma T , 

v2 

u z ■ e" = -=iTsma", (5.1) 

v2 



and 



u z x e a = —= (iiy — ia cos au x ) , 



u z x e' 



1 
1 

3 



— y= (u y — IT COS Oi' T U x ) , 



u z x e" = —= [u y +iTCOsa"u x ) . (5-2) 
V2 

Using these multiplication rules, Eqs. (|4.23|) . I|4.26|) . I|4.27|l and (|4.28l) imply, in that order, the following six relations 

^2( a 'r - a r ) = a <y - 

T 

r(a' T cos a' T + a" cos a") = (ra a cosa, 

T 

^^{K' T a' T sin a' T — K T a" sin a") — a a K a sin a , 

T 

T{e'a' T sin a' T — eq" sin a") = aa^esixia, 

T 

r (a' T e'v' T — a"ev T ) = aa„tv a , 

T 

(a' T e'v' T cos a' T + a"ev T cos a") — a a ev,jCOSot. (5-3) 



If we now use Snell's law [Eqs. (|4.16(l and (|4.17(l ]. it is easy to see that the third equation becomes identical to the 
first, and the fourth identical to the fifth. Thus, we are left with four independent equations which, for easy reference 
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in what follows, we recollect below: 

T 

^ t{cl' t cos a' T + a" cos a") = aa a cos a, 

T 

^ r (a^e'u^. - a"ev T ) = aa a ev a , 

T 

(d T e'v' T cos a!j_ + a"ev T cos a") = a a ev a cos a . (5-4) 

r 

These equations must be solved for the four unknowns a' T and a" (with r = ±), for any given incident amplitude a a . 
Before embarking on the general solution, it is instructive to consider various particular cases. 

B. Non-chiral media as a special limiting case 

If £ is zero in both sides of the boundary, so that we are referring to ordinary, non-chiral, media, the equations for 
the coefficients simplify. Neither the indices of refraction, nor the angles, depend on the sign of the polarization. In 
terms of the linear polarization amplitudes 

01 = -^=(a+ + a_), 

a 2 = -j={a+-a-), (5.5) 

and similar ones for the primed and doubly primed coefficients, the equations become 

, _ ,/ _ 1 
a i a i — ~~^ a ' 7 ' 

a 2 cosa + a 2 cosa = — -=00,0 cos a , 

i 

V2 

oieVcosa' + all ev cosa = —=oa a ev cos a . (5.6) 

Thus, the set of four equations splits into two 2x2 blocks, one for a[ and a" and the other for a' 2 and a 2 , which are 
easily solved to give 



a 2 e'v' — a 2 ev = —=<ja a ev , 



a: 



-^=a a E 1 



where 



2ev cos a 



ev cos a + e'v' cos a' 



ev cos a — eVcosa' 



2ev cos a 

E 2 = I I . (5.8) 

ev cos a' + e'v' cos a \ 11 „„„ „„„ / 1 

1 e v cosa — ewcosa 



The original coefficients for the circularly polarized amplitudes, 



a'r = -j=(a'i -i™ 2 ) (5.9) 
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and similarly for a", are then given by 

1 



2 cl I7 {E 1 +<jtE 2 ). (5.10) 



Regarding the physical interpretation, the formulas in Eq. (|5.7|) give the solution for the reflected and refracted 
amplitudes expressed in the linear polarization basis, with the corresponding quantities in the circular polarization 
basis being obtained by means of Eq. (|5.10(l . It should be noted further that, whichever basis we use to express the 
solution, we have assumed that the incident wave has a definite circular polarization specified by a, which can be ±1. 
However, we can easily obtain the solutions for an arbitrary polarization of the incident wave by superposition. If the 
incident wave is a combination of the two circularly polarized waves, then the solution is generalized to 



i^a^-Ei+oTEa). (5.H) 



or equivalently, in the linear basis, 



— / *a,o-Ei — ai E^ 
V2' 



= ^=^oa IJ E 2 =a 2 E 2 . (5.12) 

This solution in fact reproduces the well known results for the situation we are considering. It is reassuring to 
confirm, for example, that it embodies the Fresnel formulas when we take the incident wave to be linearly polarized 
along u x or u y . Let us consider the case in which the incident wave is linearly polarized along the u x direction, 
which implies that the incident wave is an equal admixture of the two circularly polarized waves, with the coefficients 
satisfying 



1 

71 

a 2 = 0. (5.13) 



a + = a_ = — ai , 



Then the solution for this case, namely 



"l 



"2 



— diEi 



= (5.14) 



in the linear basis, is immediately recognized as the so-called Transverse Electric solution. Similarly, the solution for 
linear polarization along e 2l 



a'l 



= 0. 



= a 2 E 2 (5.15) 



corresponds to the standard Transverse Magnetic solution. The Fresnel formulas are reproduced by writing the above 
results in terms of the magnetic permeability \x and the refractive index n, instead of quantities e and v. These two 
sets of quantities are related by Eq. (|3.4|) . which for non-chiral media reads 
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Eliminating e and v in favor of fi and n, the transverse electric solution can be written as 

a\ 2n cos a 



ai n cos a + cos a' 
1 M 



n cos a — cos a' 

- = 7, (5-17) 

ai n cos a + fyn' cos a' 

which, after using Snell's law to eliminate the angle of refraction, is the usual form (given, for example, by JacksonQ). 
The result for the transverse magnetic case can be reproduced similarly. 

C. Incidence from ordinary to chiral media 

A particularly simple situation, which reveals some of the features of the general solution, arises by considering the 
case in which the first medium is an ordinary dielectric, but the second one is not, i.e., 

C = 0, CVO. (5.18) 

In this case, 



n+ = n,- = n 



and by Snell's law 

1 

so that the equations for the coefficients become 



(5.19) 
(5.20) 



(a' + + a'_) - (a" + a") = a CT , 

{a' + cos a + — a_ cos a_ ) + (a" — a" ) cos a = aa a cos a , 

e'(a',v' + — a_v_) — (a" — a"_)ev = aa a ev , 

e (cl+v J- cosa , + a_v_ cosa_) + (a , + a_)evcosa = a a evcosa. (5-21) 

We can eliminate the doubly-primed coefficients in the fourth using the first, and likewise in the third using the second. 
The resulting equations are easily solved for the primed coefficients, and the doubly primed ones can be obtained by 
using the first and the second equations. Thus, the solution is easily obtained, and it can be expressed in the form 



a' i 2ev cos a 



a a \ T 



1 r , 



l -^T~) + ? £ (cos a ~ Xt cos a ' x) n ' (5 ' 22) 

where 

T r ff = (ev + Tae'v'_ T ) (cos a'_ T + to cos a) , (5.23) 

and 

t = - (T+rr - t+t-) 

= ee'v(v' + + v'_)(cos 2 a + cosa^_ cosa'_) 

+ (e 2 v 2 + e' 2 v' + v'_) cosa(cosa^ + cosa'_) . (5.24) 
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Normal incidence 



A particularly simple solution is obtained if, in addition, we assume that the incident wave travels normal to the 
plane boundary, a — 0, in which case the other angles are also all equal to by Snell's law. The above solution then 
reduces to 

- = kl + ra)F T , 
a a 2 

^ = l -{l-Ta)G_ T1 (5.25) 
a a 2 



where we have defined 

F T = 



lev 



ev + e'v' T 



€V — € V 

G T = F T — 1 = , ; . (5.26) 
ev + e'v' T 

This solution shows various interesting properties. For example, if we consider either a = + or a = — , the 
transmitted wave has the same circular polarization as the incident one, but the reflected wave has the opposite. This 
is a well known effect even in the case of ordinary media (£' = 0). The distinguishing feature in the present case is 
that the intensity of the reflected (and the transmitted) wave depends on what is the polarization of the incident one. 

If the incident wave does not have a definite circular polarization, but it is a linear combination of them, then by 
superposition the general solution for that case is 

i - 5*l>+">. 

cr 

- = JG- t £(1-t<t). (5.27) 

cr 

Let us consider the case of linear polarization along ei, as in Eq. (|5.13|1 . In that case the solution is 

< ± F 
oi ~ V2 T ' 
a" 1 

- = -j^G-r- (5.28) 
Oi V2 

Thus neither the reflected nor the transmitted waves have a definite linear polarization (e.g., a' + ^ a!_) in contrast to 
the situation considered in Section IV Bl Furthermore, while the incident wave is represented by 

lf nc \z) = ai e ie lKz , (5.29) 



the transmitted component is 



A II {z) = ^ T a 1 e' + F+e lK '+ z + e'_F^e lK '- z . (5.30) 
v2 



That is, while the incident wave is, at the boundary (z — 0), an equal admixture of the positive and negative 
circular polarizations, given by the e\ direction, the transmitted wave is split into an unequal admixture of the two 
circular polarizations. As we saw in Section IlII CI such an admixture describes an elliptically polarized wave, with 
the particular characteristics that are due to the chiral nature of the medium, as we explained there. Moreover, in 
the region z > 0, the vector that determines the direction of the electric field in the present case is obtained from the 
formula given in Eq. I|3.13|l . by making the substitutions 

A± - F±, 

6{z) l -{K'_-K' + )z, 

A(z) - \{K' + +K'_)z. (5.31) 
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D. General solution 



Other special cases can be treated similarly, such as the case of incidence from chiral to ordinary media, including 
the particular situation of normal incidence. However, we do not proceed any further along those lines, and go directly 
to give the solution to Eq. I|5.4|l for the reflected and refracted amplitudes, in the general case. For this purpose, we 
use the shorthand notation cq = cos a, c'± = cosa4 and c± = cosa± for the cosine of various angles. Then, defining 

D = {c' + c!_+clc'L){e'v' + + ('v'_)(ev++ev-) 
+ {c+c'L + c'_c'l)(e'v'_ + ev+)(e'v' + + e«_) 

+ (c'+c'i + c'_c'L){e'v' + ~ ev+)(e'v'_ - e«_) , (5.32) 

the solution is given by 

— D — (aev a +Te'v'_ T )Uv + + ev-)(ac c'_ T + Tc'lc'^) 

+ (<rev a + ev-)(e'v'_ T + T€v + )(ac c'^ + tc'_ t c'_1) 

+ ((Tev rT — ev + )(e'v'_ T — Tev~)(acoc'l_ + tc'_ t c'L) , (5.33) 

a" 

— D = (aev a + TEV- T )(e'v' + + e'v'_)(ac Q c"_ T - tc'.c') 

a a 

+ (crev a + Te'v'_ T ){e'v' T + ev_ T ){ac a c'_ T - tc' t c"_ t ) 

— (aev a — Te'v' T )(e'v'_ T — ev- T )(ac a c' T — tc'_ t c"_ t ) . (5.34) 

These formulas can be written in several alternate forms, in terms of the angle of incidence and the refractive indices, 
by using Snell's law [Eq. I|4.16|l ]. In the appropriate limits, they reduce to the special case formulas already seen. 



VI. CONCLUSIONS 



The electromagnetic properties of a medium that exhibits chirality (also called optical activity), but which is 
otherwise linear, homogeneous and isotropic, can be described in terms of the usual two parameters e and fi that 
represent the dielectric and permeability functions, and an additional parameter £ that is indicative of the chirality 
property. As we have argued here and in the references cited, such a parameterization in terms of just one additional 
function, is both complete and minimal. In the case of an infinite (unbounded) medium, an electromagnetic wave 
exhibits the phenomenon known as 'natural optical activity', which is due to the fact that the two circularly polarized 
states travel with different speed if C is non-zero. 

In this work we have considered in detail various aspects of the propagation of a wave, in a medium that is made of 
two semi-infinite media, either or both of which may be chiral, which are separated by a plane interface. Mimicking 
the procedure that is applied to the analogous problem involving ordinary dielectrics, our approach was based on 
writing the plane wave solution that satisfies the Maxwell equations in each region, supplemented by the appropriate 
boundary conditions at the interface. 

We considered first various kinematical aspects of the solution which are independent of the detailed nature of the 
boundary conditions. In particular, we obtained the Snell law and the condition for the total internal reflection effect, 
in the form that it applies to the present situation. 

In contrast with the above, the determination of the relative amplitudes of the reflected and refracted components 
is a dynamical issue that involve the details of the boundary conditions. A significant problem that we had to solve 
was precisely to elucidate what is the appropriate set of boundary conditions that must be satisfied at the interface. 
The boundary conditions depend on the Maxwell equations themselves, and since these depend on the parameter (, 
the boundary conditions necessarily involve that parameter as well. 

After finding the appropriate set of independent conditions, we solved them and obtained the amplitudes in various 
particular cases and configurations, such as one region being chiral and the other one being an ordinary dielectric, 
and the case of normal incidence. We considered various features of those solutions, and in particular we obtained in 
various cases the corresponding formula for the angle of the direction of polarization. 

In the general case, the boundary conditions form a set of four (linear) equations that must be solved for the four 
unknown amplitudes in terms of the amplitude of the incident wave. While the explicit form of the solution is not 
particularly illuminating, we obtained it and for completeness we wrote it down. 

Our work opens the way for handling a whole class of related problems involving chiral media, that we can now 
formulate in a concrete way. We can consider, for example, a slab of one material of finite thickness, inserted between 
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two semi-infinite media, and in various configurations depending on which ones are chiral or ordinary dielectrics. Pos- 
sible generalizations include also periodic or semi-periodic arrays of alternating materials, and similar arrangements. 
From an algebraic point of view, all of them will be ultimately reduced to writing down the plane wave solutions 
in each region, and imposing the boundary conditions at each interface. In this work we have shown the way for 
treating and solving all such problems systematically. In addition to the interesting and potentially important physical 
applications that we have mentioned, some of the ideas exposed here can be useful in other physics problems that 
have been studied in the literature (l8lll9|. which have a similar mathematical structure. 

Acknowledgments 

This work was partially supported (JFN) by the U.S. National Science Foundation Grant PHY-0139538. 



[1] J. D. Jackson, Classical Electrodynamics, 2nd edition (Wiley, New York, 1975). 
[2] S. Ichimaru, Statistical Plasma Physics (Addison Wesley, Massachussetts, 1992). 
[3] T. Verbiest, M. Karunanen, and A. Persoons, Phys. Rev. Lett. 82, 3601 (1999). 
[4] E. Charney, The Molecular Basis of Optical Activity (Wiley, New York, 1979). 
[5] S. Mohanty, J. F. Nieves, and P. B. Pal, Phys. Rev. D 58, 093007 (1998). 
[6] J. F. Nieves and P. B. Pal, Phys. Rev. D 39, 652 (1989). 
[7] D. R. Smith et al., Phys. Rev. Lett. 84, 4184 (2000). 
[8] E. U. Condon, Rev. Mod. Phys. 9, 432 (1937). 
[9] P. Drude, The Theory of Optics (Dover, New York, 1959). 
[10] E. J. Post, Fundamental Structure of Electromagnetics (North Holland, Amsterdam, 1962). 

[11] A. Lakhtakia, V. K. Varadan, and V. V. Varadan, Time-harmonic Electromagnetic Fields in Chiral Media (Springer, New 
York, 1986). 

[12] J. Kong, Electromagnetic Wave Theory (Wiley, New York, 1986). 
[13] J. A. Kong, Proc. IEEE 60, 1036 (1972). 

[14] C. M. Krowne, IEEE Trans. Antennas Propag. 32, 1224 (1984). 
[15] J. C. Monzon, IEEE Trans. Antennas Propag. 38, 227 (1990). 
[16] P. Hillion, Phys. Rev. E 48, 3060 (1993). 
[17] J. F. Nieves and P. B. Pal, Am. J. Phys. 62, 207 (1994). 

[18] F. W. Hehl, Y. N. Obukhov, and G. F. Rubilar, Int. J. Mod. Phys. A17, 2695 (2002), gr-qc/0203105. 
[19] A. B. Balakin and J. P. S. Lemos, Class. Quant. Grav. 19, 4897 (2002), gr-qc/0209062. 
[20] See Ref. |17j . and references therein. 



